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I .  INTRODUCTION 


With  regard  to  cases  which  do  not  give  consideration  to 
orbital  perturbations,  research  results  associated  with  optimum 
transfer  problems  between  two  circular  orbits  have  been  quite 
numerous.  Among  them,  the  most  famous  and  practically  useful  are 
dual  tangential  plans.  The  plans  in  question  are  designs  which, 
in  accordance  with  dual  impulse  methods  (dual  pulse  type),  make 
spacecraft  realize  transfer  energies  from  one  circular  orbit  (the 
initial  orbit)  to  another  circular  orbit  (the  final  orbit)  most 
economically.  The  basic  points  are  described  below.  First  of 
all,  when  the  spacecraft  moves  along  the  initial  orbit  to  the 
intersection  line  of  the  planes  of  the  two  circular  orbits  (first 
orbital  change  point),  it  enters  into  the  first  orbital  change 
(applying  the  first  impulse).  This  causes  it  to  transfer  in  o 
elliptical  orbit  (dual  tangential  transfer  orbit)  motion  with 
apogees  and  perigees  respectively  on  the  two  circular  orbits. 
After  that,  waiting  until  when  the' spacecraft  moves  along  the 
elliptical  transfer  orbit  to  the  location  of  the  terminal  orbit 
(Also  on  the  intersection  line  between  the  planes  of  “e  two 
circular  orbits.  This  second  orbital  change  point  and  the  first 

•  -t-  resoectively  set  on  the  two  sides  of  the 

orbital  change  point  are  respectively  s 

center  of  the  earth.),  it  then  enters  into  the  second  orbita 
change  (application  of  the  second  impulse),  making  it  transfer 
into  terminal  orbit  motion.  The  essential  elements  of  the  planes 
of  the  elliptical  orbits  in  question  (that  is,  angle  o 
inclination  and  ascending  nodal  point  right  ascension)  are 
deindant  on  the  key  orbital  elements  of  initial  orbits  and 
teLinal  orbits  (radius,  angle  of  Inclination,  and  ascending  node 
right  ascension).  Specific  values  are  determined  in  accordance 

with  optimized  conditions. 

However,  the  results  discussed  above  are  only 

acting  as  one  type  of  approximate  solution  ^ 

transfer  problems.  Because  the  gravitational  field 
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is  certainly  not  a  strictly  central  field  of  force,  and  space 
close  to  the  earth  is  not  an  absolute  vaccuum,  spacecraft  orbits, 
therefore,  always  exhibit  perturbations.  When  considering  key 


orbit  perturbation  factors,  optimum  transfer  problems  between  two 
circular  orbitS“-without  much  more  discussion— thereby  give  more 
practical  value  to  engineering  designs.  For  this  reason,  this 
article  draws  up— in  cases  where  flattening  effects  associated 
with  the  earth  are  considered-probes  into  optimum  plans  where 
the  greatest  energy  savings  realize-in  accordance  with  dual 
impulse  methods— transfer  from  one  circular  orbit  to  another 
circular  orbit  along  dual  tangential  elliptical  orbits. 

Obviously,  when  consideration  is  given  to  the  effects  of 
flattening  associated  with  the  earth-besides  spacecraft  moving 
in  an  eguatorial  plane  being  able  to  make  circumferential  motions 
around  the  earth-in  other  cases,  osculating  orbrts  assocrated 
with  spacecraft  motion  orbits  are  not  capable  of  marntarnrng 
round  shape.  At  most,  they  are  only  capable  of  being  near 
circular  shapes  or  osculating  orbits  which,  at 
achieve  a  round  shape.  As  a  result,  as  far  as  the 
orbits  of  this  article  are  concerned,  it  is  only  possible  to 
carry  out  solutions  by  a  type  of  approximation  in  accordance  with 
LcuLting  orbits  at  instants  of  orbital  change  being  circular 
shapes  as  well  as  with  regard  to  near  circu  ar  • 

further  rangentim  transfer  orbits 

as  a  starting  point,  on  the  one  hand,  it  is  based  on  he  orbits 

in  question  being  optimum  transfer  orbits  when  “bital 

■nrrt-  f  1  cuied  lu .  Ou  the  other  hand,  it  is 
perturbations  are  rransfer  problems 

:  ::  considerations  are  made  of  orbital  P-burbations^  es  es 
that  only  considering  flattening  effects  associated  with 
earth  certainly  does  not  mean  that  other  perturbation  factors  do 
"rexisrt  It  is  just  that,  in  order  to  simplify  discussions  of 
Torkinut  effects  or  effects  believed  to  be  created  by  other 
preturbation  factors  on  orbits,  it  is  possible  to  go  throng 
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It  is 


orbital  maintenance  systems  to  provide  their  elimination, 
also  necessary  to  put  forward  that,  with  regard  to  flattening 
perturbations  associated  with  the  earth,  consideration  is  only 
given  to  the  two  primary  effects  of  orbital  plane  forward  motion 
and  perigee  (line  of  apsides)  rotation. 

II.  DETERMINATION  OF  KEY  DUAL  TANGENTIAL  OPTIMAL  TRANSFER 
ORBITAL  PLANE  FACTORS 

If  one  wants  to  determine  optimal  transfer  plans  to  be 
realized  along  dual  tangential  elliptical  orbits  between  two 
circular  orbits  in  accordance  with  dual  impulse  methods,  one 
first  of  all  needs  to  solve  for  the  key  elements  associated  with 
the  planes  of  the  transfer  orbits  in  question. 

Assume  that  the  key  orbital  elements  of  the  two  circular 
orbits — initial  orbit  and  final  orbit  are,  respectively  as 
follows. 

Initial  orbit  radius  is  Rl.  Angle  of  inclination  is  il. 

The  ascending  node  right  ascension  associated  with  instant  t 
osculating  orbits  is  Ql^ 

Final  orbit  radius  is  R2 .  Angle  of  inclination  is  i2 .  The 
ascending  node  right  ascension  associated  with  instant  t 
osculating  orbits  is  Q2. 

Note  that  dual  tangential  elliptical  transfer  orbit 
semiparameter  is  PT.  Eccentricity  is  eT.  Angle  of  inclination 
is  iT.  The  ascending  node  right  ascension  associated  with 
instant  t  osculation  orbits  is  QT.  Perigee  angular  distance 
(measurement  from  ascending  node  to  perigee  point  along  the 

direction  of  orbital  movement)  is  a>T* 

In  accordance  with  the  meaning  of  dual  tangential,  it  is 

possible  to  know  that 
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R.+R, 


in  accordance  with  research  results  [1,2]  associated  with 
pulse  type  orbit  changes,  in  regard  to  the  two  velocity 


increments  lAF 


and 


,  they  are,  respectively. 


tAK,l=F,<^l  + 


/?,+/?,  \r,+rJ 


\avA=vA\  + 


2'-^2y-^R+R^  \r,+r 


I 

y  X  [cos/  ^cos/ ,  +  sin/  ^  sin/ ,  cos(f2  „  -  O  „ )] 


5  (3) 


'I  D  \2  /rb  M 

- ^ —  )  X  [cos/  _  cos/  2  sin/  y.  sin/  ^  cos(£J  H  )J 


rlipectively,  the'  ascending  iode  right  ascensions  associated  with 

transfer  orbits  for  the  first  orbital  change  instant  and 

respective  osculating  orbits  associated  with  initial  orbits .  nT2 
respecri  ^  ^orendina  node  right  ascensions 

and  Q22  are,  respectively,  the  ,,^hital  change 

associated  with  transfer  orbits  for  the  second  orbital  change 

instant  and  respective  osculating  orbits  associated  wit  ma 

orbits  Besides  this,  note  that  the  ascending  node  righ 

:::rc:rar:L:rirnrc::;  functloLl  relationship  between 
:guation  [3,  and  eguation  rayThlfirgr 

!isociated  with  transfer  orbit  perigee  angular  distances  do  not 
„,ake  velocity  increments  needed  for  orbital  changes  (below  we 
will  see  the  influences  on  |^F.I  values,  give  rise  to 

^^'"^''irsituations  associated  with  opting  for  the  use  of  dual 
tangential  transfer  plans  involving  the  shortest  time  periods 
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(that  is,  setting  out  from  the  first  orbital  change  point  along 
dual  tangential  transfer  orbits  and  applying  second  orbital 
changes  when  first  reaching  final  orbits),  due  to  perturbation 
effects  associated  with  movement  forward  along  the  planes  of 
relevent  orbits  given  rise  to  by  earth  flattening  and  angular 
perigee  distance  alterations,  these  two  will  lead  to  the  forms  of 
relationship  set  out  below  /15 


^22  ~  ^2l  ^  T^2 


(5) 

(6) 


In  the  equations,  and  are,  respectively,  rates  of 

change  associated  with  ascending  nodes  of  transfer  orbits  and 
final  orbits.  tT  is  the  time  gone  through  in  moving  along 
transfer  orbits  from  first  orbital  change  points  to  second 
orbital  change  points . 


(7) 


2R 


.  3.5 


cosi^/D  / s) 


(R 


1  •  -'2 
3.5 


cos/,  /  D  /  s) 


(8) 


In  the  equations.  Re  is  global  equatorial  radius;  D  is  mean 
solar  day  (using  seconds  as  measurement). 

tT  requires  considering  changes  in  perigee  (or  apogee) 
angular  distances,  that  is,  half  periods  when  orbital 
perturbations  exist  where  tT  will  deviate  from  transfer  orbits. 
Under  conditions  where  orbital  period  changes  given  rise  to  by 
flattening  earth  effects  are  ignored  as  well  as  assuming  that 
surface  areas  and  velocities  are  still  maintained  as  fixed  values 


along  orbital  motions,  tT  can  be  calculated  in  accordance  with 
the  equation  below  (in  the  vicinities  of  second  orbital  change 
points,  circular  arcs  are  taken  to  replace  transfer  orbits) 


2^  =  7’, 


l-a>r 


r: 


—  T 

1  ^  T 


180(i?,  +R,)iR,R2y 


(9) 


in  the  equation,  180  is  actually  180  degrees.  TT  is  the  period 
associated  with  transfer  orbits.  change 
associated  with  transfer  orbit  perigee  (or  apogee)  angular 

distances » 


CJ 


2R 


3-5 


Ii_=  -5 


R,  +R: 


-5cos^- )/i)  r/s) 


(1 


(10) 


(11) 


1.  •  „  i  c;  the  earth’s  gravitational  constant. 

Tn  the  equations  I  fJ-  is  tne  eai.  y 

Fron.  equation  (3)  to  equation  (11),  it  is  possible  to^^^^^ 
that  lAK.I  and  depend  on  the  key  trans 

Plane  elements  iT  and 

::^i=rtisfy  l^it  vame  equations  of  the 
forms  set  out  below 


air 


(12) 


(13) 


In  the  equations; 


f  is  the  sum  of  velocity 


with  dual  orbit  changes. 


increments 


associated 
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/=1AF,1  +  IAFj1 


(14) 


In  this  way,  it  is  possible  in 
fixed  solution  equations 


"the  end  to  obtain  the  two 

/16 


^sin/^cosj,  —  cos/j.sini,co5(ftj.,  , 

sin/^cos/ ,  -  cos/^sin,  :os(n„  -  +  sinf^sinz^sm(Q„  - j 


=  —  Sim 


s 


'^3i.  ^ 


lAF, 


V,R, 


smz 


V, 

-  sin;  ^ sin sin(Q  „  -  22 


In  the  equations; 


3t . 


R]T] 


31 


1  -  U), 


RJt 


r  360{R^+ R.KR.R.y- 


303 


i-  =  -  50^ 


IR 


31 . 


R+R. 


180(^1  +R^){R,R.) 

sin;  j. cos/ j-  /  D(\  /  s) 


30). 

(5/'"') 

31 . 


(15) 


(16) 


(17) 


(18) 


(19) 
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k  is  the  conversion  coefficient,  equaling  0.0174533  (l/°). 
This  is  due  to  angles  in  this  article  all  opting  for  the  use  of 
degrees  as  units  of  measurement.  In  equation  (18),  360  and  180 
Qj-0  really  360  degrees  and  180  degrees. 

Making  use  of  equation  (3)  to  equation  (11)  as  well  as 
equation  (15)  to  equation  (19),  through  numerical  value  solution 
methods,  It  is  possible  to  specifically  solve  for  the  values  of 
IT  and  nil  .  From  this,  specific  determinations  are  also  then 
made  for  key  transfer  orbit  plane  elements  associated  with  dual 
tangential  optimal  transfer  plans.  The  precise  determination  of 
orbit  change  points  is  discussed  below. 

III.  PRECISE  determination  OP  DUAL  TANGENTIAL  OPTIMAL  TRANSFER 
ORBITAL  CHANGE  POINTS 


After  solving  for  key  transfer  orbit  plane , elements 
associated  with  dual  tangential  optimal  transfer  plans,  use  is 
made  of  the  first  orbital  change  point,  which  should  be  on 
intersection  points  of  initial  orbit  and  transfer  orbit  planes, 
and  of  second  orbital  change  points,  which  should  be  on 
intersection  points  of  final  orbit  and  transfer  orbit  planes, 
is  not  difficult  to  precisely  determine  their  locations  [  1 . 

First  orbital  change  points,  at  angular  distances  i  .  ,  on 

osculating  orbit  planes  associated  with  initial  orbits  for  or  i 
change  instant  t  (within  the  orbital  planes  in  question,  angl 
measured  beginning  from  ascending  nodes  along  the  direction  of 
motion)  are  determined  by  the  equation  below 


sin/j.sin(nj.|  —  fij,)cos0|  + 


cos;  ^sin; , 


—  sin;  j-cosi 


;,cos(Qj., -0„)  sm0,  =0  (20) 
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Second  orbital  change  points,  at  angular  distances  162  ,  on 
osculating  orbit  planes  associated  with  final  orbits  for  orbit 
change  instant  t  (within  the  orbital  planes  in  question,  angles 
measured  beginning  from  ascending  nodes  along  the  direction  of 
motion)  satisfy  the  equation  below 

sinzVsin(Q„  -  +  [cos/^sin/^  -  sinz ^cosi^cos{Q„  -  =0  (21) 

Making  use  of  equation  (15)  and  equation  (16),  it  is  possible  to 
prove  that,  in  situations  where  orbital  perturbations  are  not 
considered,  first  orbital  change  points  are,  in  fact,  located  /17 
on  intersection  points  between  initial  orbit  and  final  orbit 
planes.  However,  when  orbital  perturbations  are  figured  in,  due 
to  the  existence  of  and  ,  first  orbital  change 

points  are  certainly  not  generally  located  right  on  the 
intersection  points  of  initial  orbit  and  final  orbit  osculating 
orbit  planes  associated  with  orbital  change  instants. 

From  equation  (20),  it  is  possible  to  know  that  has  two 
values.  From  equation  (21),  it  can  also  be  seen  that  6.  has  two 
solutions  as  well.  As  far  as  specific  selections  are  concerned, 
they  should  work  in  concert  with  in  order  to  satisfy  the 
requirements  of  dual  tangential  transfer  orbits.  With  regard  to 
this,  this  article  does  not  make  a  detailed  discussion. 

At  this  point,  dual  tangential  optimal  transfer  plans  have 

been  completely  determined. 

IV.  SPECIAL  SOLUTIONS 

With  regard  to  the  special  case  where  i2  =  0°,  from  equation 
(16),  it  is  possible  to  immediately  obtain 

(22) 
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Equation  (15)  simplifies  to  become 


1 


IR,  /'  2^, 

^  +/?7  ”  A^i  ^^^2 


(23) 


Although,  as  far  as  equation  (22)  and  equation  (23)  are 
concerned,  results  are  the  same  as  when  orbital  perturbations  are 
not  figured  in,  nT2,  however,  is  certainly  not  the  same  as  fiTl  . 

From  equation  (20)  and  equation  (22),  it  is  also  possible  to 
obtain 

sin(?,  =0  (24) 

that  is,  d  can  be  taken  as  0°  or  180  . 

From  equation  (21),  it  is  possible  to  obtain  (because  iT  /V-:- 

0° ) 


sin(0j  +Q2j)-0 


(25) 


Due  to  first  orbital  change  points  being  located  on 
equatorial  planes,  starting  out  from  limit  angles,  £221  can  be 
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taken  to  be  fill  (With  regard  to  i2  =  0°,  ascending  nodes  already 
have  no  significance.  In  all  cases,  relevent  parameters  are 
understood  from  limit  angles).  As  a  result,  from  eguation  (5) 
and  equation  ( 6 ) ,  one  has 


(26) 


Taking  equation  (26)  and  substituting  into  equation  (25),  one 
gets 

sin[0-r-(Qj.- ^2)1  =  0 

^  ^  ^  (27) 


This  is  nothing  else  than  to  say  that  9^  should  be  chosen 
from  among  the  two  values  tT  (Q^— n,)j  and  180“  tT  (Q^.  — 

in  order  to  work  in  concert  with  9^  (0“  or  180“),  making  dual 
tangential  conditions  achieve  adequacy.  It  should  be  brought  up 
that  6 2  here  is  measured  using  fi22,  calculated  in  accordance  with 
equation  (6),  as  datum.  If  one  uses  fi21  (that  is,  fill)  as 
measurement  datum,  angular  distances  9^  associated  with  /17 
second  orbital  change  points  can  provide  selection  values  as  tT  Clj 

and  180“  +  tT  flj.. 


V.  CONCLUSIONS 

Due  to  the  existence  of  orbital  perturbations,  research  or. 
optimal  transfers  between  two  circular  orbits  requires 
considering  their  influences.  When  orbital  perturbations  are  not 
figured  in,  optimal  transfer  plans  solved  for  between  two 
circular  orbits  are  only  one  type  of  approximate  solution  to 
actual  problems.  When  figuring  in  orbital  perturbations  given 
rise  to  by  flattening  associated  with  the  earth,  this  article 
gives  dual  pulse  type  optimal  transfer  plans  carried  out  along 
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dual  tangential  elliptical  orbits  (using  minimal  orbital  change 
0j20]-gi0g  act  as  an  index  of  optxmizaton)  .  Relevent  formulae 
are  capable  of  supplying  preliminary  design  uses  in  engineering. 
In  conjunction  with  this — in  regard  to  solving  dual  tangential 
optimized  rendezvous  problems  associated  with  two  spacecraft  in 
different  planes  (considering  the  effects  of  flattening 
associated  with  the  earth) — there  will  also  be  some  benefits. 
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